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A strict method to calculate the color rearrangement probability at partonic 
level in hadronic W^W~ decays is proposed. The color effective Hamiltonian He is 
constructed from invariant amplitude for the process e+e" — > W^W~ — > qiq2Q3Q4 + 
ng (n = 0, 1, 2, • • •) and is used to calculate the cross sections of various color 
singlets and the color rearrangement probability. The true meaning of the color 
rearrangement is clarified and its difference from color interference is discussed. 
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I. INTRODUCTION 

At LEP2 energy, the real W~^W~ pair production through e+e~ annihilation becomes 
possible. More accurate measurements on W mass (M^y) and other properties can be made 
at this stage of LEP project. The success of the precision measurements of Mw relies 
on accurate theoretical knowledge of the dynamics of the production and decay stages in 
e+e^ W^W~ 5'i^2?3^4- However, the possible Color Rearrangement (CR) may obscure 
the separate identities of two W bosons so that the final hadronic state may no longer be 



considered as a superposition of two separate W decays. Thus the W mass determination 
may be distorted by the color reconnection effect in hadronic W^-pair decays. This effect was 
first studied by Gustafson, Pettersson and Zerwas (GPZ) It attracts a lot of studies in 
recent years Considering that the two Ws decay into two quark pairs qiq2 and ^3^4, 

GPZ assume that these two original color singlets can be rearranged into two new ones ^1^4 
and ^3^2 with a probability (|), and then energetic gluons are emitted independently within 
each new singlet, which implies that the CR occurs before the parton shower process begins. 
But Sjostrand and Khoze do not regard this instantaneous scenario as a very likely one 
The reason is that the decay vertices of two W bosons are in general separated in space-time, 
and therefore the hard gluons (with Eg > T^) are produced incoherently by the two pairs 
^1^2 and ^3^4 0,3. So there are two color singlets Ci and C2, each containing a qq pair and 
a set of gluons. Furthermore they argue that the CR in Perturbative QCD (PQCD) phase 
only comes from the color interference which should be very small. Hence they conclude that 
the non-perturbative contribution dominates the CR effect because the two color singlets Ci 
and C2 coexist later during the relatively larger space-time scale of hadronization compared 
with that of W^'s life. The non-perturbative CR probability is controlled by the space- 
time overlaps of the color field induced by two groups of partons in Ci and C2. Later on, 
Gustafson and Hakkinen stress that CR can only originate from the partonic level and argue 
that the hard gluon emission will enlarge the CR effect. Because there are increasing ways of 
color recombination between the partons of Ci and those of C2 with the growing number of 
emitted gluons, even though the rearrangement probability is only for each way, the total 
probability at the partonic level may be greatly enhanced to the order ~ (/ + l)(m + 1)/N^ 
where / and m are the numbers of gluons in Ci and C2 respectively. The final probability can 
in principle be modified by multiplying the factor of order ~ (/ + l)(m + 1)/N^ by unknown 
functions of variables which characterize the space-time overlaps of the color fields induced 
by partons of Ci and C2. But Gustafson and Hakkinen's analysis on the total probability at 
the partonic level is only a qualitative one. It does not include many other ways of forming 
singlets and the calculation is not based on a strict formulation. 
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We should keep in mind that the color fields stretched between partons must be restricted 
in a preconfined state, i.e. a singlet. Thus in PQCD phase, CR means the transformation 
from a set of original singlets to that of new ones. Hence the CR probability from Ci (con- 
taining 5152) and C2 (containing qsq^) to new recombined ones where qiq^^ and ^3^2 belong 
to different singlets should and can be estimated by PQCD with more reason and accuracy. 
In this paper, we try to calculate the rearrangement probability from a strict systematic 
approach of PQCD. This approach is based on the color effective Hamiltonian He which 
is built from the recursive form jsHlOl of the invariant amplitude M for the process: 

e+e" ly+iy^ qiq2mi + ng, n = 0, 1, 2, ■ ■ ■ . (1) 

The color effective Hamiltonian He is used to calculate the cross sections and the probabilities 
of various rearranged color singlets formed by final partons. The physical significance of our 
approach lies in that it includes all of effects caused by the different space-time intervals 
between the decay vertices of two VT-bosons. The CR probability we obtain shows that the 
CR of PQCD stage is not negligible. This seems different from what Sjostrand and Khoze 
conclude in ref. 0,0]. In our approach, the meaning of CR can be clearly defined and the 
difference from color interference can be easily elucidated. 

The outline of this paper is as follows: in section Jl|, we give the invariant amplitude 
M„ for the process (0) in the recursive form; then the color effective Hamiltonian He is 
abstracted from M„ in section |IIT| thirdly, we use He to analyze the color singlet structure 
of the parton states q'i^2?3^4 + ^9 with n = 0, 1, 2 and give the rearrangement probabilities 
in section |^ finally, a summary is given in the last section. 

II. CROSS SECTION 

The differential cross section dcxn for the process (|^) is 

dan = $|M„prfp„+4(gi, Q2, Q3, Q4, ^1, ■ ■ ■ , Kn), (2) 
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where $ = ^ is the flux factor multiphed by a spin average factor, and the phase space 
factor dpn+4, is deflned by 

here P+, P_, Qj (i = 1, ■ ■ • , 4) and -ft^j (j = 1, ■ ■ ■ , n) denote the 4-momenta of e+, e~, 
quarks and gluons respectively. 

According to Feynman rules, the matrix element M„ of the process ([l|) (see flg.l) can be 
written as 

m=0 V=Ue,-y* ,Z0 ^3"^ 
y<Si,r{Qi; Ki, ■ ■ ■ , Kjn', Q2)S^i{Qz'-, Km+l, ■ ■ ■ , Kn', Q4), 

where denotes the 4-momenta of boson, and D'^'^{W^'^) is the propagator of 
boson; £^ is the polarization tensor of leptons; 5',^ is the current containing quarks and 
gluons which depends on the 4-momenta, helicities and color indices of the outgoing partons. 
In the recursive form, the current S,^ can be expressed by 

S,{Q,; K,, ■ ■ ■ , K^; Q2) = legT E (^"'^"' " " " 1, 2, ■ ■ ■ , ^; Q2), (4) 

P{l,2,---,m) 

where Qs is QCD coupling constant; = ^ and A" is Gell-Mann matrix for S'f/c(3), and 
5*,^ is the spinor current (for detail, see refs. 
Substituting eq. (|D into eq. (^, we obtain 

n 

m=0 P(l,...,m) P(m+1, ■■■,»!) 

where the indices i, k (j, /) denote the color (anticolor) of outgoing quark (antiquark), and 



y-Syi{Ql] Ki, ■ ■ ■ , Km', Q2)S^i{Qz] Km+l, ■ ■ ■ , -ft'n! Qi)- 



can be written as follows 



(6) 

=T;(P+)[^err]t.(P-)Z^(^)"^(g'^)[^eFX^(Q',W^+,W^-)], V = 7*,^°, 



where the repetition of indices represents summing (we use this convention unless exphcitly 
specified); is the fermion-boson vertices for the boson 5 = 7*, Z", W^; Q = — P+ 
and Q' = -P+ + P-', D^°'^{Q''^) is the propagator of the vector boson V (=7* or Z^)] 
F^ui^iQ'^ W') is the three-boson vertex of YW+W- {V = 7*) or Z^W+W- {V = Z^). 

Up to now, we use the recursive formula to write down the invariant amplitude M„ for 
the process e~^e~ W^W'^ — > 919293^4 + "^Q- this form, we see that M„ can be clearly 
expressed as a uniform color part multiplied by a momentum function of partons. The 
effective Hamiltonian can be found from this form of the amplitude. This is what we 
shall do in the next section. 



III. COLOR EFFECTIVE HAMILTONIAN He 

In ref. 0, from PQCD, a strict formulation has been proposed to calculate the cross 
section of color singlets for the process e+e^ ^ jZ^ qq + ng. Now we use the same 
approach to abstract the color effective Hamiltonian H^. for the process (|I]) from the invariant 
amplitude M„ given in eq. (|^). Then He is found in the following form: 

n 

-'-'c V-' ^ jjy-^ Ji^^(qigi---gmq2)(q3gm+i---9uq4) 

m=0 P{l,---,m) P{m+l,---,n) 

(7) 

m=OP(l,...,m)P(m+l,-,n) 

^"^(l?igi---9m92)(<?39m + l ■■■97194)' 

where = {R^ , B^)u and \E'{+ = {R^ , B^)u are the color creation operator for 
quark and antiquark, respectively. The color octet operator Gu of gluon u is defined by 

Gj,^. = -^(A-A-)} = v[/^+vI/+ - ivl;^+vl;+ 55, (8) 

where A"" (a^ = 1, ■ ■ ■ , 8) is color operator for gluon u and is defined in ref. 0. Here He is 
another expression of S matrix, so it is not necessarily Hermitian. 
For a final color state \f >, its cross section can be calculated by 
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ai = f\H,\Q > pdp„+4, (9) 

where |0 > is the initial state free of color. If 

|/>= |[vl/+vi/i+(Gi...Gj;.][x[/+vl/!,+(G^^^ >, 
after summing over all of the color indices, we have 

E4 =I^J:\< f\Hc\0 > ?dpn+A = / $ < mtHc\^ > 

f f (10) 

= /$|M„pdp„+4 = an- 

It shows that the calculation of the ordinary cross section via He returns to the original 
form. So the validity of for the process (|I]) is verified. 

For the parton system qiq2<lz<ii + ^9-, the color state is composed of the color charges of 
(111 5 5'3 5 ^4 ^'^d n gluons. It belongs to the color space 

3i(g)3;(g)33(g)3:(g)8i(g)---(g)8„. 

There are many ways of reducing this color space. Corresponding to each reducible way is 
one set of orthogonal singlet sub-spaces whose bases contribute a complete and orthogonal 
set of color singlets. For a color singlet set >, /c = 1, ■ ■ ■}, we have 

|/fc></fc| = l, <fk\fi>=Ski (11) 

and 

EcT^ =I^E\< fkmo > \^dp^^^ = / $ < 0\H+\fk >< fk\H,\0 > dp^+, 
k k (12) 

= / $ < 0\H+He\0 > dpn+A = I $|M„|2rfp„+4 = an. 

This is the result of unitarity. 

According to eq. (|^), for instance, one can get the concrete expressions of color effective 
Hamiltonian Hen for the the process (|l]) with n = 0, 1, 2 as follows: 
a.) for the process e~^e~ W^W~ '?i^2Q'3^4 = 0), 

Heo = (vl/+vl/^+)(vl/+v[;^>)X(,,,^)(^3-^); (13) 
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b. ) for the process e~^e — > W^W — qiq2<l'i<lAgi ('^ = 1); 
where 

c. ) for the process e+e^ W^W~ (1iQ2Q3'1a9i92 = 2), 

Hc2 = + + ifc2) 

where 

Hc2 = \ E (*i^2^)(^3jG'iA:'^l*4^)-^(gig2)('?39l9294)' 
P(l,2) 

P(l,2) 

IV. COLOR SINGLET STRUCTURE OF THE FINAL PARTON SYSTEM 

In this section, we try to use He derived in last section to study the color singlet structure 
of the final parton system, give the CR probability and compare our results with those of 
other authors. As the gluon number grows, it becomes more and more difficult to calculate 
(T„, and the number of different ways of forming color singlets increase drastically. To 
illustrate what our approach is and how it works, we study only three lowest order cases: 

(1\1Mai Q1Q2Q3Q491 and qiq2q3q49i92- 

A.) For the parton system q'ig2Q'3^4) due to the reasons given in ref. we have no 
need to consider the color configurations resulting from the reductions of 3 3 = 3* 6 
and 3* 3* = 3 © 6*. Hence the color space 3i 32 Sa (8) 34 has only two reducible ways 
which correspond to the color configurations with and without CR. In the following, we will 
discuss these two cases in detail. 
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(14) 



(15) 



For the reducible way 

(3l (g) 3*2) (g)(33 (g) ^ (I12 812) (g)(l34 834) ^ (I12 (g) I34) 0(8l2 (g) 834) 

the color singlet set is {|/o >, i = I, 2} where 



1/0 >= Imi^l^mj'^t) >, 1/0 >= ^\Tr{G^,Gu) >, (16) 



where G'^y^ = ^ti^y'^ - {xy = 12, or 34) denotes the color octet state formed by 

Qx and Qy. In |/q >, '^'^^'^^2^ and '^^j'^'!^ represent the two initial color singlets within ^1^2 
and ^354. Defining the probabilities of |/q > (z = 1, 2) as = I^\<fo\yo>\^'ip^ ^^^^^ 



is the cross section for e+e^ W^W~ C[i(i2Q3Q4 Hco is given in eq. (|T3D, we find that 
Pq = 100% and Pq = 0. This set corresponds to the original color configurations. 
For the other reducible way 

(3l (g) 3,1) (g)(33 (g) 3;) ^ (li4 814) (g)(l32 832) ^ (ll4 (g) I32) 0(8l4 (g) 832), 

the color singlet set is {|/q >, i = I, 2} where 



1/0 >= TTl(^i,^4+)(^3^,^r) >, 1/0' >= —\Tr{G^,G,,) > . (17) 



Their probabilities are given by 

1 !<^\<f^\H,o\0>\'dp, 1 2 /$! </o^|gco|0> pc^P4 8 

(To y (To 9 

One sees that the probability of |/q >, where qiq^ and 53^2 form two singlets \l/5''j\l/4"^ and 
^37-^2^5 is |. This state is just the CR case first discussed by GPZ [|l|. But here it clearly 
shows that the states with and without CR belong to two different completeness sets. In 
the new set, the state /q is the color singlet made up of two color octets which are formed 
from color charges of the pair 51^4 and ^3^2 respectively. So all of the four partons must in 
principle hadronize as a unity in this state, though we do not know how to treat this kind 
of hadronization rigorously. Note that the state |/o > is different from |/d >. But in LUND 
string fragmentation picture, in |/g > which is built up by two octets, the neutral color flow 
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connects qi with and with q^, and the substrings stretched within the pair gigg and 
^3^4 are treated as subsinglets (in fact, at the partonic level, they are not singlets, but only 
color neutral objects). Then the hadronization result of |/q > certainly has no difference 
with those of |/q > in LUND model. 

B.) The color space 3i 83 {8> Ss 84 (g) 81 of the parton system gig2Q'3^45'i can be reduced 
in two ways: one corresponds to the case with no CR, and the other to that with CR. Here 
we discuss the latter case. Corresponding to the following reducible way 

(3i 3^) {8)(33 (8) 3^) (8) 81 ^ (ii4 e 814) 0(i32 e 832) (8) 81 

^ (li4 (8) 832 (8) 81) e(8i4 (8) I32 (8) 81) e(8i4 832 (8) 81), 
the completeness set of color singlets is {\fl >, j = 1, • • • , 4} where 



1/1 > = ^|Tr(GMGi)(vi/+ >, 1/2 > = _^|Tr(G32Gi)(vI/+ >, 
Iff > = yX|Tr({GM,Gi}G32) >, 1/1' > = ^|Tr([GM,Gi]G32) >, 



(19) 



with 



{Gi, 02)1 — G\iG2k + G\iG\f. - lTr{GiG2)Sl, 
Their probabilities {i = 1, ■ ■ ■ , 4) are given by 

Pi = < n\H,i\o > \'dp, = < fi\Hi,\o > + < n\H'^,\o > i^dp, 

= ^/ $[E I < n\HL\0 > p + 2Re{< 0\H',^\n >< /{|^ci|0 >)]dp5, 



(20) 



where cxi is the cross section of the process e~^e — > W^W — > qiq 213(1 4,91, and Hd is given 

in eq. (|T3|). They are related to the energy and PQCD parameters (e.g. and 
2 

etc.). / $ S I < /il-f^ci|0 > pc?p5 are proportional to ai. Thus in the last line of eq. (|20|), the 
first term does not depend on these quantities. The second term is the color rearrangement 
caused by the color interference. Even if it vanishes, the color rearrangement still exists. 
This shows that the color interference contributes to only part of CR. Additionally, it is easy 



to verify J2 Pi = ^ which is the natural result of eq. (|T2D 

1=1 
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We notice that in \fl >, Tr^Gi^Gi) is the color singlet formed by gi, ^4 and gi, while 
^3x^2''' the color singlet by ^3 and gg- The case of \ ff > is similar to that of f^. These two 
states are just the CR states discussed by Gustafson and Hakkinen in their naive model [Q. 
But the meaning of the probabilities of these two states in their work is different from that 
in ours. For the process e+e~ W^W^ '?i^293^4fl'i) they consider two possible cases: in 
one, Qi is radiated from gigj with a probability P, and in the other, gi is radiated from ^3^4 
with a probability (1 — P). There are no interference terms. For each case, they give the CR 
probability |. So according to their analysis, one can derive that the total CR probability 
is I for the whole process, since |P + |(1 — P) = |. But in our calculation, the probability 
is expressed in eq. (^O]). One cannot distinguish which source the gluon is radiated from. 
Since in eq. (0), the first term H^^ describes the case that the gluon is emitted from 
and the second one H^^ describes that the gluon is from W~ . These two terms combined 
to give probabilities in eq. (|^). Note that in this paper, as was done in ref. 0], we also 
only consider the hard gluon emission. In this case, we find that the interference terms are 
negligibly small, and we obtain Pi = P^^ ~ |, Pf ^ ^ and Pf ~ | from eq. (|20|). We notice 
that \fl > 01 \ ff > represents the color singlet formed by two color octets: one is the octet 
formed by two octets from gi and qiq^ which are symmetric in \fl > or antisymmetric in 
1/^ >, and the other is from 5352- Iii this two states, the four quarks are included in a whole 
singlet so they do not emerge in two different subsinglets. As we recall that CR means the 
two original singlets Ci (containing ^1^2) and C2 (containing ^3^4) are rearranged to make 
qiq^ and ^3^2 enter into two new different subsinglets, our result of CR probability is 2/9. 
It is interesting to further look at how LUND model treats the hadronization of the states 
Iff > and 1/1 >. According to the neutral color flow picture of LUND model, such states 
result in two neutral color flows: one from qi to ^2 via gi, and the other from ^3 to q^, or one 
from qi to ^2 ? a-^id the other from ^3 to ^4 via gi . These two neutral flows were approximated 
to singlet string pieces which fragment into hadrons independently. So in LUND model, no 
difference exists of the hadronization result of |/f > and \ ff > from those of the states with 



10 



no CR. In this sense, the CR probabihty is also about |. 

C.) For the parton state (?i^2'?3^45'ifl'2, to obtain the singlet set that results in the CR, 
we reduce its color space 3i ® 83 (8) 84 (g) 81 82 as follows: 

(3i (8) 31) 0(33 3*2) 81 82 ^ (li4 e 814) 0(l32 e 832) 81 82 

^ [ii4 132 81 82] e[ii4 832 81 82] e[8i4 132 81 82] (21) 
e[8i4 0832 081082]. 

Obviously the color singlets corresponding to the first three terms in the above equation 
lead to CR. The last term of eq. ( pTD can be further reduced as, e.g. 

(a) . (814 081)0(832 082) ^ 08. ® 183.08 J e [08^,08, 00^3. 08 J' 

(b) . (8i4082)0(83208i)^ [l8^^(g,8^0l83^(g,8je[O8^^(g,8^0O;^^(g,3j, 

(c) . (814 832) 0(81 82) - [ig,, (g, 18, (g, 8.] e [08,, 83. (2) o;^ 3 J , 

(d) . 814 0(832 081082) ^ 1814 0(832 81 082), etc., 

where 08^^081 (^832 82 ^ "denotes the nonsinglet state formed by qi, and gi {q^, ^2 and 
(72), and 08^4 081 ^^832 082 denotes the color singlet formed by this two nonsinglet states, 
and so on. The reduction ways (a) and (b) lead to the same color configurations. The 
cases (c) and (d) give a slightly different total CR probability from that of the case a or b, 
because all of the singlets reduced in (c) and (d) are not color rearranged ones, while the 
first singlets reduced in (a) and (b) contribute to the total CR probability. As we can see 
in the following, the difference of total CR probability between this two groups of reduction 
cases is about 8%. Hence the CR probability slightly depends on the reduction way one 
chooses. But our current knowledge of QCD is not enough for us to determine which the 
physical reduction way or the physical singlet set is. As an example, here we discuss the 
reduction way (a). The corresponding color singlet set {I/2 >, j = 1, ■ ■ ■ , 7} is given by 

1/2 >= ^|(^i.^D(*3,^DTr(GiG2) >, I/I >= ^|(vl/i,vl/^)Tr(G32{Gi,G2}) >, 

I/I >= ^|(*i,n)Tr(G32[Gi,G2]) >, I/I >= ^|(vl/3.vl/^2)^KGi4{Gi,G2}) >, ^^^^ 

I/I >= M'^3^^)TriGu[GuG2]) >, |/| >= ||Tr(Gi4Gi)rr(G32G2) >, 

\f >=^l(O8..08.®o;^,08.)>> 
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where is the normahzation constant for |/J >• The probabihties of I/2 > = 1, ■ ■ ■ , 7) 
are defined by 

^^ Jt|<m.|0>l'dp. ^ , = 1, ...,7, (23) 

0-2 

where 02 is the cross section for the process e+e" W^W^ (I1Q2Q3Q49192, and Hc2 is 
given in eq. (|T5l). Note that Hc2 depends on momentum configurations of partons, so a2 and 
P2 are functions of ^/S and the PQCD parameters. In fig. 2, we give the probabihties Pi at 
y/S = 170 GeV. One can notice that they are not sensitive to Ymin- 

From eq. (|22D , we see that in I/2 >, '^u'^l, '^sj'ifi and Tr{GiG2) are three color subs- 
inglets within 51^4, ^3^2 and gig2 respectively. Obviously this color separated singlet as we 
call it is not covered by the model of Gustafson and Hakkinen. Its probability is only about 
1.3%. In 1/2^ > and |/| >, ^u^i is the color singlet within ^1^4, while Tr ((732 {Gi, 6*2}) 
and Tr{G32[Gi, G2]) are the singlets formed by three octets: one is symmetric and the other 
antisymmetric for gi and g2- The situations of |/| > and |/| > are similar to those of 
I/I > and I/I >. In |/| >, Tr^GuGi) {Tr{G-i2G2)) is the singlet formed by gi, and 
Qi (?3) Q2 and (72)- According to our definition, at \fS = 170 GeV , the total CR probability 
for e+e^ — > W^W~ Q'i^25'3^45'i5'2 is about 28% (see fig. 2) for the singlet sets (a) and 
(b), and about 20% for (c) and (d). Following a similar discussion as we have made for 
e'^e^ W^W~ (li(l2^'i^A9i^ see that the meaning of the probabilities of these states 
are different from what Gustafson and Hakkinen imply in their estimation of the CR prob- 
abilities. From their analysis, the total probability of color rearrangement for this process 
is ~ I at the partonic level. 

Note that Gustafson and Hakkinen regard the CR probability at the partonic level as 
the upper limit which occurs when the two decay vertices of and W~ coincide in space- 
time. The final value should be decreased by taking into account that two W^-bosons decay 
at different space-time points. But in our approach, the CR probabilities at the partonic 
level already include the effects of non-overlapping decay vertices of two VT-bosons. This 
might be the main difference between these two approaches for the meaning of the CR 
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probability at the partonic level. 



V. SUMMARY 

The study of CR in hadronic W~^W~ decays is of significance for both precisely measuring 
the mass of W and clarifying the vacuum structure of QCD. The main goal of this paper 
is to provide a strict approach to deriving how many CR singlet sets exist in each final 
parton system and calculating the CR probability at the partonic level. To meet this goal, 
we use the recursive approach to give the invariant amplitude M„ for the processes e'^e~ — > 
— > gig2Q'3^4 + f^9 {n — 0, 1, 2) first; Then from M^, we abstract the corresponding 
color effective Hamiltonian He] Finally, He is applied to calculate the CR probability. We 
find that the CR probabilities are § and about 20% ~ 28% (at ^/S = 170 GeV) for 
e"^e~ — > W^W~ — > qiq^Q^QA + ''^9 with n = 0, 1 and 2 respectively. 

Summarily, the following points should be noted: 

• Our result of the CR probability in PQCD stage is an accurate one at the tree level. It 
already contains the effect caused by different space-time intervals between the decay 
vertex of W'^ and that of W~ because our approach is a matrix element method. It 
shows that the CR probabihty in PQCD phase is not quite small. Our approach to 
studying CR seems different from Sjostrand and Khoze's. The difference may lie in the 
different definition of CR at the partonic level. We defined CR as the transformation 
from the original color singlet set to a new one where qiq^ and 93^2 belong to different 
color subsinglets, while Sjostrand and Khoze defined that as the "antennae" 14 and 
32 in momentum space. The "antennae" ij — ^ ^'^'^ where pi, pj and k are the 
momenta of qi, qj and the gluon radiated from the dipole qiqj. 

In our approach, the color Hamiltonian H^ is abstracted by factorizing the invariant 
amplitude M into the uniform color part multiplying the momentum one. Then we 
choose and explore one color singlet set which contains the CR singlets where qiq^ and 
?3?2 belong to different color subsinglets. By doing projection of CR states to H^ i.e., 
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< fcR\Hc\0 > where fcR denotes the CR singlet, we derive the CR amphtude, then the 
cross section and the probabihty. In the procedure of calculating the CR probability, we 
see that even the color interference terms vanish, there are also CR contributions which 
are not of interference origin, which shows that the color interference contributes to 
only part of CR. In Sjostrand and Khoze's approach, all color indices are summed over 
in |Mp; the "antenna" terms of (14) and (32) imply CR only arise in the interference 
sector of |M|^. The "antennae" can be considered as the sources of further dipole 
cascading , a parallel description of parton cascading. As we remember, the parton 
shower process can be well described by the CLAP equation where the color indices 
are summed over at each branching point, so CLAP equation is a probability evolution 
equation. Another parallel scenario is based on BFKL evolution equation |jl3|] where 
the color configuration is also smeared. BFKL equation is closely related to the color 
dipole model. One great advantage of these approaches compared with our matrix 
element method is that they are easy to be implemented by Monte-Carlo simulation. 
But in these approaches, the other CR sources besides interference cannot be fully 
revealed. 

We are not sure whether the final event properties predicted from our approach 
is different from those from Sjostrand and Khoze's, because in order to produce final 
hadrons from the partonic states one has to apply a phenomenological model to de- 
scribe the hadronization process of such color states as |/o >, > and |/J > etc.. 
Rigorously say, such states are beyond the scope of the currently used fragmentation 
model. 

• Our results for three lowest order cases correspond to three certain reduction ways of 
the color space which lead to one specific color singlet set for each case that contains 
the CR states. For q'ig2?3^4 5'i^29'3^4fl'i case, the corresponding total CR probability 
that we obtained is the same as what Gustafson and Hakkinen estimated. This is 
a curious coincidence, because for each of these two cases, only one complete and 
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orthogonal singlet set leads to CR, i.e., all of the CR states belong to the same singlet 
set and orthogonal to each other. For q'ig293945'i5'2 case, there are four singlet sets 
that lead to CR. For two of them the total CR probabilities are the same (~ 28% 
at y/S — 170 GeV), while for the other two, they are a httle smaller (~ 20% at 
\fS = 170 GeV). The singlet set we choose is one of the former. It implies that in 
general, for a certain final parton system, there may exist different CR sets and the 
corresponding CR probabilities. We are far from making the physics choice among 
them. On the other hand, in this case, some of the CR states are not orthogonal to 
each other since they do not belong to the same singlet set. For example, one of the 
rearranged singlets, |/' >= \\Tr{Gx4G2)Tr{G^2Gi) >, belongs to a singlet set that 
is different from the set of |/| >, so they are not orthogonal to each other, as we 
see < /'I/I >= I 7^ 0. The probability of a certain CR state is defined in its own 
set and is meaningless outside it. But Gustafson and Hakkinen do not discriminate 
different color singlet sets in their estimation where the probabilities of non-orthogonal 
singlets are added up together and the total CR probability is more than necessary. 
However, our result displays the same tendency that the gluon emission enlarges the 
CR probability as predicted by Gustafson and Hakkinen. 

• The CR effects observed in the final hadron events still depend on the hadronization 

mechanism and the vacuum structure, but they originate from the color reconnection 
at the partonic level. Only when the CR properties at the partonic level are made 
clear is it possible to investigate the nature of the hadronization mechanism and the 
vacuum structure. 



15 



REFERENCES 

[1] G. Gustafson, U. Pettersson and P. M. Zerwas, Phys. Lett. B 209 (1988) 90. 

[2] T. Sjostrand and V. A. Khoze, Z. Phys. C 62 (1994) 281. 

[3] T. Sjostrand and V. A. Khoze, Phys. Rev. Lett. 72 (1994) 28. 

[4] G. Gustafson and J. Hakkinen, Z. Phys. C 64 (1994) 659. 

[5] L. Lonnblad, Z. Phys. C 70 (1996) 107. 

[6] Yu. L Dokshitzer, V. A. Khoze, L. H. Orr and W. J. Stirling, Nucl. Phys. B403 (1993) 
65. 

[7] Qun Wang and Qu-Bing Xie, Phys. Rev. D 52 (1995) 1469; QCD Argument on Color 
Flow Model, Qun Wang, Qu-Bing Xie and Zong-Guo Si, to be pubhshed in Phys. Lett. 
B. 

[8] F. A. Berends and W. T. Giele, Nucl. Phys. B 306 (1988) 759. 
[9] F. A. Berends, W. T. Giele and H. Kuijf, Nucl. Phys. B 321 (1989) 39. 
[10] F. A. Berends and W. T. Giele, Nucl. Phys. B 294 (1987) 700. 

[11] Physics at LEP2, CERN 96-01, Vol.1, Editors:G.AltareUi, T. Sjostrand and F. Zwirner, 
Pages 124 and 190. 

[12] G. Gustafson, Phys. Lett. B175 (1986) 453; Nucl. Phys. B392 (1993) 251; B. Andersson, 
P. Dahlqvist and G. Gustafson, Phys. Lett. B214 (1988) 604. 

[13] Ya. Ya. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822; E. A. Kuraev, 
L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45 (1977) 199. 

Figure Captions 

Fig.l: e+e" W'^W" Q'i^293^4 + ng {n = 0, 1, ■ ■ ■) process. 
Fig. 2: The probabilities Pi as a function of Ymin- 
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